Supercooled Liquids under Shear: Computational Approach 
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Abstract. We have performed molecular dynamics simulations for a model two-dimensional soft-core mixture in a super- 
cooled state. The mixture exhibits a slow structural relaxation in a quiescent state, however, the relaxation is much enhanced in 
sheared states. There observed surprisingly small anisotropy both in the coherent and incoherent density correlation functions 
even under extremely strong shear which is 10 3 times faster than the structural relaxation rate. The present simulation results 
agree well with predictions of the recently developed mode-coupling theory in shear. 



As liquids are cooled toward the glass transition, the 
dynamics is drastically slowed down, while only small 
changes can be detected in the static properties. One 
of the main targets in theoretical investigations on the 
glass transition is to identify the mechanism of the drastic 
slowing-down. Beside this fundamental problem, a strik- 
ing example occurs when one brings glassy materials 
away from equilibrium, for instance, by changing tem- 
perature rapidly or applying shear flow to them. There 
appears a variety of unique phenomena such as aging 
or shear thinning. Although these phenomena are not 
only conceptually but also practically important, physi- 
cal properties of glassy materials in nonequilibrium con- 
ditions has not yet been understood well. A couple of 
years ago, we performed extensive molecular dynamics 
(MD) simulations in two dimensions (2D) and three di- 
mensions (3D) on binary soft core mixtures with and 
without shear flow. It was found that the dynamical prop- 
erties of the mixtures under shear can be fairly mapped 
onto those at quiescent states at higher temperatures 
We found also the surprising isotropy in the tagged 
particle motions even under extremely high shear, which 
may justify the simple mapping idea. In the present 
study, we calculate intermediate scattering functions by 
using the method proposed by Onuki [3,4] to investigate 
microscopic dynamics of glassy materials in shear flow. 
Simulations have been done in 2D to compare the present 
computational results directry with the theory developed 
recently in 2D ||5[ . 

Our model system is composed of two different parti- 
cle spices 1 and 2, which interact via the soft-core poten- 
tial 

v a p(r) = £(<J a p/r) n , O a p = {<J a + Op)/2, (1) 
where r is the distance between two particles, and a, j3 S 



1,2. We take the mass ratio mi/mx = 2, the size ratio 
O2/01 = 1-4, and the numbers of particles N\ = N2 = 
5000. Simulations were performed in the absence and 
presence of shear flow with being particle density and 
temperature fixed at n = [N\ +N2)/V = 0.8/of and 
kgT = 0.526e, respectively. In the sheared case, we kept 
the temperature at a constant using the Gaussian con- 
straint thermostat to eliminate the viscous heating effect. 
Our method of applying shear is as follows: The system 
was at rest for t < for a very long equilibration time 
and was then sheared for t > 0. Here we added the aver- 
age velocity jy to the velocities of all the particles in the 
x direction at t = and afterwards maintained the shear 
flow by using the Lee-Edwards boundary condition. Sim- 
ulation data have been taken and accumulated in steady 
states which can be realized after transient states. 

Figure 1 (a) shows the geometry of shear flow in the 
present simulation. As shown in Figure 1 (b), shear flow 
with the rate 7 advect a positional vector r as 

r(t) = r + ytr y e x (2) 

in the time duration t, where e a is a unit vector in a G x,y 
axis. A similar advection can be defined in Fourier space 
for a wave vector k as 

k(t)=k-ytk x e y . (3) 

The above definition enable us to calculate Fourier com- 
ponent k of the time correlation function 

C(k,t) = (A^(0)A k{t) (t)) (4) 

in shear flow. We calculated the incoherent and coherent 
scattering functions by using the definitions 

Fs{Kt) = ^£ e [-«W)-r.-W-k- ri (0)}]^ (5) 
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FIGURE 1. (a) Geometry of shear flow, (b) Shear advection 
in real space, (c) Shear advection in Fourier space, (d) Sampled 
wave vectors (t=0). 



and 

F ( k) r) = 1/ £<.[-*«*»(')] £ e [*Ty(0)]\ (6) 

respectively. To examine anisotropy in the dynamics, the 
wave vector is taken in four different directions kio, koi, 
ki i , and k_n, where 

k ai3 = ^4==(ae* + j3e y ), (7) 

fco = 2%/ai, and a,/3 6 0, 1 as shown in Figure 1 (d). 

Figure 2 (a) and (b) show Fsi(k.,t) and Fn(k,f), re- 
spectively, at r = 0.526 with and without shear flow. 
Here the subscript 1 denotes the smaller particle com- 
ponent. The so-called a relaxation time T a of the present 
mixture is defined by 

F 11 (k ,'Ca)^Fs 1 (ko,'C a )=e- 1 (8) 

in the quiescent state. The followings have been found; 
i) Fs\(k,t) and F\i(k,t) behave quite similarly both in 
quiescent and sheared states, ii) Shear accelerates drasti- 
cally the microscopic structural relaxation in the super- 
cooled state. The structural relaxation time i a decreases 
strongly with increasing shear rate as x a ~ 7~ v with 
V ~ 1 . iii) The acceleration in the dynamics due to shear 
occurs almost isotropically. There observed surprisingly 
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FIGURE 2. Intermediate scattering functions under shear; 
(a) incoherent part and (b) coherent part. 



small anisotropy in the correlation functions even in ex- 
tremely strong shear jx a ~ 10 3 . This simplicity in the 
dynamics is quite different from behaviors of other com- 
plex fluids such as critical fluids or polymers in shear. 
Finally we note that the recent mode-coupling theory in 
shear flow [5] predicts almost the same behaviors. 
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